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ABSTRACT 


Initially  nonconducting  gas  is  ionized  by  a  thin  viscous 
shock  wave.  Upstream  there  can  be  no  magnetohydrodynamic  inter¬ 
action  because  of  the  zero  conduct ivi ty,  but  the  conducting 
downstream  region  may  have  a  magnetic  structure  which  interacts 
with  the  flow  variables.  A  theoretical  analysis  is  made  in  the 
zero-magnet  ic-Prandtl  -number  (,!non- viscous" )  limit,  i.e..  Ohmic 
dissipation  is  the  dominant  diffusion  mechanism.  Unlike  magne¬ 
tohydrodynamic  shocks  in  a  pre- ionized  gas,  ionizing  shock 
waves  are  not  necessarily  plane-polarized.  Thus  "skew"  shock 
structures  can  exist,  in  which  the  upstream  3nd  downstream 
magnetic  field  vectors  and  the  shock  wave  normal  do  not  all  lie 
in  a  single  plane.  The  existence  of  the  viscous  subshock  at 
the  front  of  the  magnetic  layer  requires  the  upstream  Mach 
number  to  be  greater  than  unity.  Explicit  solutions  are  given 
for  typical  values  of  the  governing  parameters,  showing  how  the 
magnetic  field  vector  rotates  about  the  shock  wave  normal  as 
its  transverse  component  changes  in  magnitude  through  the  shock 
layer.  Skew  shocks  are  necessarily  sub-Al fvenic  downstream. 
Unlike  the  pre- ionized  case,  the  range  of  trans-Al fvenic  shock 
waves  is  not  excluded,  since  these  shocks  can  absorb  Alfven 
waves  within  their  structure.  With  strong  magnetic  fields  it 
is  possible  to  achieve  very  high  downstream  temperatures  by 
Joule  heating.  Alte-natively,  in  some  cases,  magnetic  energy 
can  be  fed  into  directed  kinetic  energy,  producing  an  overall 
expansion  shock. 
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JL*  Introduction 

• 

This  article  concerns  the  magnetic  structure  of  ionizing 
,  shock  waves  when  Ohmic  dissipation  is  the  dominant  diffusion 

mechanism.  Viscous  and  heat  conduction  effects  are  neglected 
in  the  magnetic  shock  layer,  and  are  important  only  in  an  in¬ 
finitely  thin  subshock,  which  is  assumed  to  exist  at  the  front 
of  the  magnetic  layer.  This  subshock  is  responsible  for  ioniz¬ 
ing  the  previously  non- ionized  upstream  gas  so  that  it  becomes 
an  electrical  conductor  somewhere  within  the  subshock,  and  re¬ 
mains  conducting  everywhere  downstream.  As  will  be  shown  in 
§  2.3,  the  detailed  variation  of  the  (scalar)  electrical  con¬ 

ductivity,  although  important  in  the  full  analysis  of  the  prcb- 
lem,  can  be  suppressed  by  a  suitable  nonuniform  "stretching" 
t  of  the  normal  space  coordinate.  The  initial  analysis  is  thus 

equivalent  to  that  of  a  hypothetical  gas  which  has  zero  conduc¬ 
tivity  upstream  of  the  subshock,  and  constant  (finite)  conducti¬ 
vity  everywhere  downstream.  The  neglect  of  subshock  fine-struc¬ 
ture  means  that  the  present  analysis  is  the  "cuter  expansion"  in 
a  matched  asymptotic  expansion  procedure  as  described,  for  exam¬ 
ple,  by  Van  Dyke  (1964).  The  inner  expansion  and  matched  compo¬ 
site  solution,  together  with  the  analysis  of  the  stretching  ef¬ 
fects  of  variable  conductivity  (due  to  variations  in  the  degree 
of  ionization  and  in  the  electron  temperature)  will  be  presented 
at  a  later  date. 

*  Several  parameters  govern  the  behaviour  of  ionizing  shock 

waves.  Naturally  it  is  not  feasible  to  present  details  for  a 


full  range  of  all  the  parameters.  For  this  reason,  only  cases 
of  infinite  upstream  Mach  number  Mj  wi 1 1  be  considered,  i.e., 

[flow  kinetic  energy!  »  [thermal  energy],  upstream.  For  Mach 
number*-  above  Mj  ~  20,  this  is  certainly  a  good  approximation, 
because  the  thermal  terms  which  are  neglected  in  the  analysis 
are  O(Hj^). 

The  shock  wave  normal  defines  one  direction  (e.g.,  the 
x  axis,  say).  Perpendicular  to  this,  the  constant  transverse 
electric  field  defines  another  suitable  reference  direction 
(e.g.,  the  z  axis,  say).  The  other  mutually  perpendicular 
transverse  direction  (y  axis)  completes  a  Cartesian  coordinate 
system.  The  special  cases  of  "normal"  shock  waves  (no  upstream 
transverse  magnetic  field  components)  and  "oblique"  shock  waves 
(the  upstream  magnetic  field  has  finite  x  and  y  components, 
but  no  z  component)  will  be  considered  in  a  subsequent  paper 
[Leonard,  (;970b)].  "Skew"  shock  waves  (in  which  the  upstream 
magnetic  field  involves  ai 1  three  components)  will  be  discussed 
here,  but,  primarily,  for  only  one  set  of  upstream  angles  (see 
§  3).  The  term  "skew"  seems  to  be  appropriate  for  this  latter 
type  of  shock:  whereas,  in  normal  and  oblique  shock  waves,  the 
upstream  and  downstream  magnetic  field  vectors  and  the  shock 
wave  norma!  lie  in  a  single  plane,  in  skew  shocks  they  do  not. 
Hence  skew  shocks  represent  ' r,e  most  general  type  of  geometrical 
conf igurat ion.  The  possibility  of  such  shock  waves  was  mentioned, 
in  passing,  by  Cowley  (1967),  who,  however,  confined  his  atten¬ 
tion  to  the  more  commonly  studied  cases  involving  only  a  single 
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transverse  component.  As  will  be  shown  in  the  companion  paper, 
in  oblique  shocks,  the  magnetic  field  vector  lies  in  a  single 
plane  throughout  the  entire  shock  structure,  provided  the  con¬ 
ductivity  is  a  scalar.  Hence  oblique  shocks  are  always  plane- 
polarized.  This  is  also  true  of  normal  shocks.  The  structure 
of  “transverse"  ionizing  shock  W3ves  (no  normal  magnetic  field 
component)  is  discussed  in  Appendix  1  of  that  paper,  where  it 
is  shown  that  such  shocks  must  always  be  of  the  magnetica! ly- 
trivial  gas  shock  variety,  when  Ohmic  dissipation  is  the  only 
significant  diffusion  mechanism. 

it  will  be  shown  in  §  3.3  that  skew  shocks  are  necessarily 
sub-A I fvenic  downs trean.  Skew  shocks  which  are  entirely  sub- 
Alfvenic  always  have  a  range  of  structures,  uniquely  determined 
by  the  strength  of  the  transverse  electric  field,  which  is  taken 
to  be  a  free  parameter.  Unique  structures  can  also  be  found  for 
a  range  of  ( t rans- A1 fven i c )  shock  waves  in  which  the  ups t ream 
region  is  super-Al fven ?c,  i.e.,  un i ty- local -Aifven-number  occurs 
within  the  shock  structure.  The  fact  that  these  t rans- A1 fven ic 
ionizing  shock  waves  are  stable  to  rotational  Alfven  disturbances 
(in  constrast  to  the  anomalous  trans-Al fvenic  "solutions"  in  the 
p re- ion ? zed  case)  has  been  discussed  previously  [Leonard,  (1970a)]. 
In  terms  of  typical  ionizing  shock  structure  behaviour,  there  is, 
in  fact,  no  significant  distinction  between  the  trans-Al fvenic 
and  completely  sub-Al fvenic  regimes. 

In  the  next  section,  the  governing  equations  will  be  devel¬ 
oped.  The  behaviour  is  governed  t?y  two  first-order  differential 


equations  for  the  transverse  magnetic  field  components,  with 
auxiliary  variables  being  defined  algebraically.  An  integral 
coordinate  transformation  decouples  the  basic  equations  from 
variable-conductivity  effects.  In  §  3,  typical  behaviour  of 
the  solutions  Is  discussed  with  particular  reference  to  the 
parametric  variation  of  the  electric  field  strength.  Optimi¬ 
zation  of  the  electric  field  for  the  achievement  of  maximum 
downstream  temperature  is  also  discussed.  it  appears  that  the 
conversion  of  magnetic  energy  into  thermal  energy  by  Joule  heat¬ 
ing  can  be  responsible  for  raising  the  downstream  temperature 
by  a  factor  of  several  times  the  downstream  subshock  temperature. 
Under  certain  ci rcumstances,  magnetic  energy  can  be  converted 
into  directed  kinetic  energy  in  the  magnetic  layer,  so  that  after 
an  initial  (viscous)  compression,  a  re-expansion  takes  place 
downstream,  giving  the  possibility  (for  relatively  large  upstream 
magnetic  pressure  ratios)  of  overall  expansion  shocks. 

In  the  examples  of  §  3  and  the  Appendix,  a'l  results  are 
obtained  by  exact  numerical  evaluation  of  algebraic  expressions  or 
by  numerical  integration  of  the  governing  differential  equations. 
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2.  Governing  equat ions 

2. 1  Phys ical  1 aws 


It  is  assumed  that  a  shock-fixed  coordinate  system  can  be 
chosen  such  that  the  shock  structure  becomes  steady  and  one¬ 
dimensional.  Oimensional  variables  carry  a  bar.  The  x 
coordinate  is  taken  as  the  single  independent  va  iable,  which 
is  positive  downstream,  with  its  origin  at  the  subshock.  Uni¬ 
form  conditions  (d/dx  -*  0)  are  assumed  far  upstream  (x  =  -») 
where,  for  the  moment,  completely  general  conditions  are  assumed 
for  the  velocity  Vj  =  [vx],v  ,  ,vzj],  magnetic  field  = 

C^l  »Byj  »B2i  3,  and  electric  field  =  £EX1  »^yl  »^zl *  A  tran¬ 
sition  to  uniform  downstream  conditions  (as  x  -*  +<=>)  given  by 
¥2*  anc*  ¥2*  's  sou9^t*  The  most  symmetrical  form  of  the 
governing  equations  seems  to  arise  if  the  coordinate  system  is 
chosen  such  that  the  downstream  ¥  field  and  velocity  are  in 
the  same  plane  [the  (x,7)  plane,  say].  This  can  always  be 
done  as  follows: 


(i)  Rotation  about  the  x  axis  -*  B2  =  ^BX2»®y2,0^; 

(ii)  Translation  parallel  to  the  z  axis  -  y2  =  -''x2,^y2’^^' 
since  a  Galilean  transformation  does  not  change  the  magnetic 
field  in  a  nonrel at i vi st ic  theory. 


A  translation  parallel  to  the  y  axis  is  still  available. 
This  is  chosen  so  _hat  the  upstream  velocity  can  be  taken  as 
v,  =  [v  ,,  0,  v  ,]  for  convenience.  In  general  the  upstream 
magnetic  field  may  have  three  nonzero  components.  The  form  of 
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the  electric  field  can  be  determined  from  Ohm's  law  for  the 
conduction  current  density  ]J, 

J"  =  a[F  +  v  x  F],  ( 1  ) 


where  a  is  the  scalar  electrical  conductivity.  Ampere's 
equation  in  the  steady  state  gives 


x  (B/u)  =  1  +  pcv  ~  "J, 


(2) 


where  u  is  the  magnetic  permeability,  and  pc  is  the  negli¬ 
gibly  small  charge  density.  Applied  in  the  downstream  uniform 
region,  where  it  is  assumed  that  "c? 2  £  0,  equations  (1)  and 
(2)  require 


F«  —  Va  x  Ba  , 


(3) 


so,  in  the  chosen  coordinate  system,  F2  =  [0,0,Fz23.  Farav'ay's 
law  in  the  steady  case  requires  vxF  =  0;  thus  the  transverse 
components  of  IT  are  constant  throughout  the  shock  layer: 

Ty(x)  3  0,  and  Fz(x)  3  Fz2=  ^z»  say.  At:  this  stage,  no  preli¬ 
minary  information  is  available  on  F  (Poisson's  equation 
merely  implies  that  the  charge  density  will  be  extremely  small), 
hence  the  electric  field  can  be  written  as 


F(x)  =  CFx(x),  0,  Fz3. 


(4) 
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The  solenoidal  condition,  v*F  =  0,  requires  the  constancy 
of  the  normal  magnetic  field  component,  F  (x)  3  F  =  const. 

The  transverse  components  are  variable,  hence 

B(x)  =  [Bx,  Fy(x),  F2(x)].  (5) 

Viscous  and  heat  conduction  terms  will  be  neglected  unless 
the  gradients  of  velocity  and  temperature  become  '‘infinite11 
(i.e.,  when  these  variables  change  by  order  unity  in  a  length 
which  is  very  small  compared  with  the  "magnetic"  length  scale). 
This  is  the  zero-magnet ic-Prandtl -number  approximation: 

Pm  =  tiov  -0,  (6) 

where  v  is  the  kinematic  viscosity.  In  this  limit,  the  Ohmic 

diffusion  length  resulting  from  the  nonzero  electrical  resistivity 
_  1 

(=  o~‘  ),  is  larger  than  diffusion  lengths  of  the  relatively 
small  viscosity  or  heat  conductivity.  Mathematically  this  leads 
to  the  previously  mentioned  outer  expansion  in  a  matched  asymp¬ 
totic  expansion  procedure  in  terms  of  a  small  *  parameter  Pm. 

It  must  be  stressed  that  the  discontinuities  appearing  in  the 

*  In  the  author's  experience  in  other  analyses  (Leonard,  1969  ), 
matched  f 1 rst-order  composite  solutions  give  very  "reasonable" 
results  even  when  the  perturbation  parameter  is  as  large  as 
one-thi rd. 


?gw- 
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present  solutions  can  (in  principle)  be  resolved  by  using  suit¬ 
able  inner  expansions  and  matching  to  form  uniformly  valid 
composite  solutions;  but  this  wi  11  not  be  attempted  here. 

For  steady  flow,  the  mass,  momentum,  and  energy  conserva¬ 
tion  equations  can  be  written,  respectively 

I*  Tpv )  =  0,  (7) 

v.fpyv  -  BB/u  +  [p  +  Ef2/  (2u)]«3?  }  =  0,  (8) 

__  _  r\ 

V. [(pu  +  ^Pv  )v  +  pv  +  (rxl)/u]  =  0,  (9) 

where  p  is  the  density,  p  the  pressure,  and  u  the  specific 
internal  energy.  Electrostatic  terms  have  been  ignored  for 
reasons  explained  by  Sherd  iff  (1965).  For  convenience,  ideal 
gas  state  and  internal  energy  relations  are  assumed,  introducing 
the  effective  gas  temperature  T  : 

p  =  pRT, 

and 

u  =  RT/(y-1), 

where  R  is  the  gas  constant  and  y  the  specific  heat  ratio. 


As  used  here,  T  -  Th  +  aTg,  where  "h"  refers  to  heavy  particles, 
"e"  to  electrons,  and  a  is  the  degree  of  ionization.  Note  that 
this  is  not  the  overall  temperature,  T  =  T/(l+a). 

03  * 


(10) 

(ii) 


]] 


ionization  energy^  Has  been  neglected  because  only  large-Mach-number 
cases  will  be  studied. 

2.2  Nondimens ional  equations 

Appendix  1  of  a  subsequent  paper  [Leonard,  (1970b)]  considers 
the  case  of  purely  transverse  ionizing  shock  waves  fB  =  0). 

Thus  it  is  necessary  here  to  consider  only  those  shocks  for 
which  "B  4  0.  This  component  is  therefore  used  in  the  nondimen- 
s  ional  izat  ion  of  the  following:  B  =  B/B  ,  v  =  v/' v  , ,  E  =  £/  (v  ) 

***  A  ***  A  I  **  _Xarl  A  I 

P  =  ”p/ p | ,  and  T  =  T/'f , .  Equations  (7),  (8),  and  (9)  take  the 
integrated  one-dimensional  forms  [ using  equations  (10)  and  (11)]. 

(12) 

(13 

(14) 

(15) 


pvx  =  const. [lj  *  1, 

(this  has  been  used  in  the  following) 

v  +  — —7  +  (B^  +  B^)  =  const. [1 3 

x  vxY„2  V 

-  1  +  ^+  +  Bzl>- 

2 

vy  '  axlBy  =  const.Cn 
v  -  rtx1Bz  =  const. [2]  =  0, 


2 

"•This  term,  ~a/M],  could  be  formally  kept,  but  the  (weak)  coupling 
with  a  introduces  mathematical  complications  which  disappear  when 
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Mf(Y-l) 


+  ?(v2  +  v2  +  V2)  -  a2,E  B 
- '  x  y  z  xl  z  y 


=  const. [1 3 


~T  +  +  axl  BzlJ  “  axlEzByl»  (16) 


Hf(y-l) 


where  the  x  component  magnetic  pressure  ratio  is 


“xl  -  (5f/M)/(»,v2,)  =  I/Af, 


(  Aj  being  the  shock  Alfven  number),  and  the  shock  Mach 
number  Mj  is  given  by 


v"?  =  v^/(RT,) 


Note  that  equation  (15)  requires  the  z  components  of  the 

velocity  and  magnetic  field  to  be  proportional  throughout  the 

shock  layer.  In  particular  this  must  be  true  upstream;  this 

2  4  2 

fact  has  been  used  in  replacing  v2j  by  ax]B2j  on  the  right 
hand  side  of  equation  (16).  The  constants  appearing  in  these 
fluid  equations  have  been  evaluated  in  the  regions  shown  in 
square  brackets;  they  are  seen  to  depend  on  the  parameters: 
y,  M| ,  axj,  B  j,  Bzj,  and  Ez.  Equations  (13)  —  (16)  are  in¬ 
validated  if  gradients  of  velocity  and  temperature  become  very 
large  as  in  the  viscous  subshock,  however  these  equations  must 
be  satisfied  immediately  in  front  of  and  behind  the  subshock. 


I* 
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If  Is  eliminated  between  equations  (1)  and  (2)  the 
nondlmens lonal  component  equations  become 

Jx  =  Rm(Ex  +  vyBz  -  vzV  *  -  pcvx  “  °»  (,9) 

Jy  =  -  dBz/dx  =  Rm(v2  -  vxBz),  (20) 

and 

J2  =  dBy/d*  =  ^^2  +  VXBy  “  Vy)»  (2*  ) 

where  J  =  J/ [TSx/ ( iiIR )  ] ,  x  =  x/TR,  and  a  quasi- local  magnetic 
magnet  ?c  Reynolds  number  Rm(x)  =  [uvxj"£r]c(x)  has  been  intro¬ 
duced  using  a  reference  length  which  is  as  yet  unspecified. 

[  Rm(x)  can  be  regarded  as  a  nondimens ional i zed  conductivity.] 

In  equation  (19)  the  extremely  small  charge  density 
p  =  p  /[I  /  (uZ  v  ,)]  has  been  included  for  later  reference, 

C  C  X  K  X 1 

but  has  been  completely  neglected  in  equations  (20)  and  (21). 

Recall  that  ~o  (and  hence  Rm  )  is  proportional  to  the 
degree  of  ionization  and  inversely  proportional  to  the  electron- 
heavy  particle  collision  frequency, which  depends  on  the  local 
(non-equilibrium)  electron  temperature.  Tins  coupling  with 
ionization  and  electron  temperature  effects  (which  may  occur 
over  a  length  scale  comparable  to  the  magnetic  length  because 
of  ionization  lags  and  non-equilibrium  Joule  heating  of  the 
electrons)  can  be  suppressed  by  the  coordinate  transformation 


of  §  2.3.  Hence  it  is  not  necessary  at  this  stage  to  give 
explicit  equations  describing  the  electron  temperature  behaviour 
and  the  ionization  dynamics.  Note  however  that  it  wi 1 1  be  as^_ 
signed  that  a,  which  is  zero  upstream,  does  have  some  finite 
value  everywhere  on  the  downs tream  side  (x  >  0+)  of  the  vis¬ 
cous  subshock.  !n  addition,  precursor  ionization,  if  it  exists, 
will  be  assumed  to  be  confined  to  a  length  scale  comparable 
to  the  viscous  length,  or  else  to  be  of  such  small  magnitude 
that  significant  magnetic  effects  do  not  occur  upstream  of  the 
subshock.  (The  upstream  change  in  magnetic  field  is  proportion¬ 
al  to  the  precursor  ionization  length  scale  t imes  some  average 
degree  of  ionization;  thus  the  solutions  of  the  present  analysis 
are  applicable  only  to  those  cases  for  which  this  product  is 
smal 1 . ) 

Substituting  equations  (14)  and  (15)  into  equation  (19) 

gives 

By|  )b2  -  ax!BzBy^  “  c-  (22) 

i  .e. , 

Rm[Ex  -  (axlByl)Bz]  *  0.  (23) 

Upstream  (where  Rm  =  0)  this  is  automatically  satisfied. 

Note  that  equation  (20)  requires  Bz=  const. =  B2j  everywhere 
upstream,  but  equation  (23)  still  allows  Ex(x<0)  to  be 


arbitrary.  Downstream,  Rm  is  nonzero,  so  equation  (23)  requires 


Ex  “  <^I8yl)B2-  ;>0-  (2,,) 

and  from  Poisson's  equation  the  charge  density  is  propor- 

v* 

tional  to  the  derivative  of  this,  which  in  turn  is  proportional 
to  Jy  [from  equation  (20)3.  To  an  excellent  approximation 
the  nondimens ional  relationship  is 

pc  =  (*xl/£,2dEx/d*  "  (Vxl/p)2°xlBylJy’  (25) 

where  c  is  the  electromagnetic  propagation  speed.  This  extremely 

small  charge  density  is  caused  by  insignificant  deviation  of 

£  from  the  value  given  by  equation  (24)  [solve  for  p  in 
x  c 

equation  (19)3. 

In  the  general  case,  Bj  =  [1,B  |,B2|3,  V|  =  [1 ,0,ay|Bz|3, 
equations  (20)  and  (21)  require  that  upstream  (where  Rm  =  0) 
and  across  the  subshock:  B(x^0+)  s  Bj .  Since  IB  does  not  vary 
upstream,  there  can  be  no  magnetohydrodynamic  interaction  on  the 
fluid  variables,  so  v(x<0  )  =  v,;  but  note  that  across  the  sub- 
shock,  x  =  0  -*  0+,  v  changes  discontinuously.  In  the  down¬ 
stream  magnetic  layer  B  and  v  may  change  continuously  in 
magnitude  and  in  both  direction  angles,  finally  approaching  the 
uniform  downstream  state  in  which  B  j  =  v  2  =  Clearly 
neither  the  magnetic  field  vector  nor  the  velocity  vector  lies 
in  a  single  plane  in  the  general  case  of  skew  shocks. 
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The  above  discussion  assumes  that  the  electrical  conducti¬ 
vity  is  a  scalar.  When  Hall  currents  are  present  in  conducting- 
upstream  shock  waves,  as  described  by  Leonard  (1966),  £  and  v 
contain  spiralling  components  in  the  (y,z)  plane,  although 
their  upstream  and  downstream  limits  must  lie  in  a  single  plane 
with  the  shock  wave  normal.  This  spiralling  behaviour  is  also 
expected  to  occur  in  ionizing  shocks  (whether  skew,  obi ique  or 
normal),  when  Hall  currents  are  important,  but  will  not  be 
considered  here,  i.e.,  it  will  be  assumed  that  a  is  a  scalar. 

Equations  (20)  ano  (21)  are  the  fundamental  "magnetic11 
equations  for  the  structure  of  the  transverse  magnetic  field 
components  in  the  downstream  shock  layer.  The  coupling  with 
the  transverse  velocity  components  can  be  eliminated  explicitly 
by  using  equations  (14)  and  (15),  giving 


-ZT  =  ~  (axl-  vx)By-1’ 

dx  xi  x  y 


(26) 


where 


7T  =  -  vx)Bz’ 


E~  -  Ez  +  °*lByl  - 


const. 


(28) 


It  is  notationally  convenient  to  leave  the  coupling  with  v 

X 

in  these  equations  although,  in  fact,  v  =  v  (B  ,B_,}  in  the 

A  A  y 

magnetic  layer,  which  means  that  equations  (26)  and  (27)  are 
implicitly  coupled. 

If  Vy,  vz,  and  T  are  eliminated  from  equations  (13)  — 

(16),  the  relationship  between  v  ,  B  ,  and  B  can  be  found; 

x  y  z 

in  an  implicit  form  it  is 


f(vx,By,Bz) 


„2  rvvx 

ax 1 _y- 1 


*xl][<VByl>(VByr2V  +b2 


Bz.3 


+  7^17  <V  1  Hvx-  V*)  =  o. 


where 


bv4  =  Byt(vx)  =  ta2,-  YVx/(y-l)]  ’ 


and  where  the  downstream  subshock  velocity  v"  is 


v*  =  vx(x=0+)  =  (y-1  +  2/M| )/  (y+1  ) . 


Note  that,  everywhere  upstream  and  across  the  subshock, 

B  =  B  |,  and  8=6,,  hence  equation  (29)  requires  v  =  1 
Y  y  *  z  z  *  x 

(for  x  *  OJ,  or  vx=  v*  (for  x  =  0+). 
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Equation  (29)  represents  a  surface  f  =  0  in  (v  ,B  ,B  ) 

x  y  z 

space,  !t  can  be  rearranged  as 


(By  -  Byt)2  +  B2 


2  2  *Y+1  *  ^vx-1  *  ^vx~vP 

(By,—  B  4)2  +  B2,  +  — rf  rrf  * 

axl  taxl (Y-1  H'vx:i 


In  this  form,  it  can  be  seen  that  planes  vx=  const,  intersect 

the  surface  in  circles  with  centres  at  B  =  B  i(v  ).  B  =0. 

y  yf  x  '  z  * 

and  of  radius  equal  to  the  square  root  of  the  right  hand  side 
of  the  equation  (32).  From  equation  (30)  the  "centreline" 
curve  is  seen  to  be  a  hyperbola  in  the  (v  , B  )  plane  with  an 

a  2  x  y 

asymptote  at  vx  =  axi  (y— O/Y.  The  radius  also  has  a  singularity 
at  this  value.  Note  that  f  =  0  contains  the  upstream  point 
Cvx=l,  By=Byj,  ancj  the  subshock  downstream  point 

tVX=Vi-  W  Bz=Bzl]- 

In  §  2.3  a  new  space  coordinate  x  is  introduced  which 
eliminates  the  variable-Rm  factors  in  equations  (26)  and  (27); 
thus  By(x)  and  Bz(x)  can  be  found  by  direct  simultaneous 
(numerical)  integration,  the  normal  velocity  vx(x)  being 
computed  from  equation  (29)  as  the  integration  proceeds.  The 
density  is  p(x)  =  l/vx,  from  equation  (12).  Equations  (14) 
and  (15)  give  the  transverse  velocity  components  v  (x)  = 
ax1  (By  —  By.j ) ,  and  vz  =  axiBz*  Equation  (13)  can  be  rearranged 
to  give  the  temperature 
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T  =  vxfl  +  vMf[l  -  vx  -  (33) 

2  2  2 

where  Br  =  +  B2.  A  decrease  in  the  transverse  magnitude 

Bp  corresponds  to  a  conversion  of  magnetic  energy  into  thermal 
energy.  The  nondimens ional  pressure  and  internal  energy  are 
given,  respectively,  by 


P  a  P/(Plvxj)  =  t/(vxyM?)’  (34) 

and 

u  s  u/v^}=  T/ CyM^(y-I)].  (35) 


For  steady  flow,  and  in  the  absence  of  viscous  dissipation 
and  thermal  conduction  (i.e.,  in  the  magnetic  layer,  but  not 
in  the  subshock),  the  dimensional  entropy  equation  given  by 
Sherd  iff  (1965)  becomes 


(36) 


where  all  the  entropy  production  is  due  to  Ohmic  dissipation. 
The  nondimens ional  form  of  this  equation  is  discussed  in  the 
next  section  where  the  entropy  s  is  given  as  an  explicit 
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integral  involving  B^,  Bz>  vx,  and  T, 

Finally,  it  is  sometimes  convenient  to  rescale  the  temper¬ 
ature  in  terms  of  the  viscous  temperature  T*  =  T(x=0+),  which, 
from  equation  (33),  with  B  =  B  ,,  and  v  =  v*  is 

\  i  «  y\ 

T*  =  v*[l  +  vMf(l-vx)] 


=  2v(v-1) 

(Y+l)2 


HfO- 


2yM 


(y — 1  )M| 


,1 


Then  the  rescaled  temperature  is 


T  *  T/T* 

(y+D2vxC1  -  vx  +  1/ (YMf )  -  ia^TB^-B^)] 

2 { Y — 1  ) C 1  -  (y-1)/(2yM^)]{1  +  2/ [ (v— 1  )M2 ] 3  ' 


(37) 


(38) 


a 

In  this  form,  the  Mach  number  dependence  of  T  is  0(1). 


2.3  Integral  coordinate  transformation 

This  analysis  assumes  that  the  electrical  conductivity 
a  is  identically  zero  from  x  =  -  ®  to  x  =  0_  (just  upstream 
of  the  viscous  subshock)  and  that,  from  x  =  0^  (just  down¬ 
stream)  to  x  =  +“,  a  is  finite  and  will,  in  general,  vary  with 
x  (depending  on  the  local  values  of  th.-  degree  of  ionization 
a  and  the  electron  temperature  ).  Nondimens ional 1 y : 
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Rm  =  Rfn(a,Te),  where  a  and  Tg  are  functions  of  the  space 
coordinate.  Formally,  equations  (26)  and  (27)  can  be  rewritten 


uu  A 

3x  ’  E*  ~  (,xr  vx>V 


"  9y(VV* 


-  (a",-  vx)Bz 


*  92<v82)'  say> 


where  a  new  space  coordinate  x  has  been  introduced  satisfying 
the  differential  relation 


ax  =  Rm  dx,  for  x,  x  >  0+, 


locally.  Equations  (39)  and  (40)  can  be  integrated  [using  the 
auxiliary  equation  (29)]  to  give  8^  =  B^(x)  and  Bz  =  Bz(x), 
independent! y  of  the  variation  in  a  and  T  .  The  analysis  is 
completed  by  finding  a  and  Te  as  functions  of  the  new  coord i 
nate  x  (using  appropriate  equations,  not  discussed  here); 
then  equation  (41)  is  integrated  (numerically)  to  give 
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=  x  (x )  =  fX 


which  can  then  be  inverted  to  give  x  =  x(x).  When  this  func¬ 
tion  is  used  in  the  arguments  of  the  dependent  variables,  the 
latter  are  thus  given  as  functions  of  the  original  nond imens ional 
space  coordinate  x;  e.g.,  By  =  By[x(x)],  etc. 

The  only  other  equation  to  be  affected  by  the  coordinate 
transformation  in  the  present  case  is  that  for  the  entropy, 
equation  (36).  This  equation  can  be  nond imens i ona 1 ? zed  to  give 


^  *  gW 


2  2  Rm(9?  +  9?> 

— f— =- 


where  s  =  s/R,  and  and  gz  are  given  by  equations  (39) 

and  (40).  The  coordinate  transformation  dx  =  Rm~^dx  may  now 
be  applied  and  the  Rm  -  factors  cancelled,  giving,  on  integra¬ 


tion. 


2  2 

0  0  pX  (g  +  g  ) 

=  s(x)  =  yMjcxx j j  — y-  T  2  dx  +  s(0+). 


Note  that  the  factor,  =  (B^/p)/ (pj RTj )  =  gives  a 

measure  of  the  ratio  of  the  upstream  magnetic  pressure  to  the 


The  gas  subshock  entropy  jump  gives  s(0+)  =  0(log  Mj). 
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gas  pressure.  However,  for  large  Mach  numbers,  the  downstream 

temperature,  being  the  same  order  of  magnitude  as  the  viscous 

o  2 

temperature  I',  is  0(Mj),  so  that  s  ~  axj  rather  than 
.  In  this  case  it  is  convenient  to  use  the  rescaled 
temperature  T,  giving  [from  equations  (37)  and  (38)],  for 

AsMAG  5  5  “ 

A 

7  2 

For  small  values  of  axl»^sMAG  =  — ^°x1 ^ ’  ^ut  ^or  *ar9e  va^ues 
7*2 

ax|,  T  ~  axj  [from  equation  (38)],  so  that  As^g  ~  j)(l)  in  this 
case. 

2.4  Mathematical  formulation 

The  problem  to  be  studied  in  this  paper  is  summarized  here, 

a 

for  reference.  Shock  structure  profiles  of  v,  B,  T,  and  As^g 
are  required  as  functions  of  the  transformed  space  coordinate  x. 
No  attempt  will  be  made  to  find  x(x)  in  this  report. 

Given  explicit  values  of  the  governing  upstream  parameters 
Cy,  Mp  axl,  Byj ,  Bzl,  Ez],  the  values  of  v*  and  E"  can  be 
computed  using  equations  (28)  and  (31).  The  basic  problem  is  to 
compute  the  downstream  (x  a  0+)  magnetic  structure  B(x), 
using  equations  (39)  and  (40)  [and  the  fact  that  Bx  =  1  ]  with 
equations  (29)  and  (30)  as  auxiliary  equations  for  vx(By,Bz). 

The  initial  conditions  are  those  just  downs tream  of  the  subshock 


«M<>c<9v  *  9z)/f]dx 

2(V-1)[1  -  (y-1)/(2yM^)H1  +  2/Cfy-  1  )MfD} 


(x  =  0+):  vx  =  v*t  By  =  Byl ,  B2  =  Bzl.  Except  in  isolated 
special  cases,  (numerical)  integration  of  the  two  coupled  first- 
order  nonlinear  differential  equations  for  By  and  Bz  can 
proceed  from  the  subshock  towards  downstream  infinity,  at  which 
point  the  system  (usually)  has  a  nodal  singularity.  The  other 
dependent  variables,  v(x),  T(x),  p(x),  p(x),  u(x),  and  s(x), 
can  be  computed  from  the  respective  algebraic  or  integral  relations 
given  by  equations  (29),  (33),  (12),  (34),  (35),  and  (44),  or 
the  equivalent  relations  for  the  thermodynamic  variables  involv- 

A 

ing  the  rescaled  temperature  T. 

Note  that  the  downstream  singular  point  [  (d/dx)+oo~*  0], 
which  must  1 ie  on  f  =  0,  is  found  from  equations  (39)  and  (40) 
to  be  given  by  the  intersection  of  this  surface  with  the  hyper¬ 
bola  g2=  gy=  0  [i.e.,  with  Bz2=  0,  By2=  E*/(axl  -vx2)]. 

Since  this  point  is  in  the  (v  ,B  )  plane,  it  is  found  by  the 

x  y 

intersection  of  the  hyperbola 


(46) 


with  the  plane  curve  f  =  0  given  by 

fc(VV  =  f<VBy0) 

=  °xl  [^T  •  °xl  ][  (By'“yl  )  (By+Byr2By4)  '  Bzl  3 

+  T$T  (vx  -  1  )  (vx  -  vi>  “  °>  <47> 


frcxn  equation  (29). 
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3.  Skew  shock  structure 


3 . 1  Typical  behaviour 

In  (vx*By*Bz)  space,  the  structure  trajectories  of  skew 

shocks  start  at  (1,B  jf3z|),  which  is  the  upstream  sing'  *ar 

point  corresponding  to  x  =  This  point  is  on  the  surface 

f  =  0.  Note  that  dB^/dx  =  dBz/dx  =  0  at  this  point  (and, 

in  fact,  everywhere  upstream,  x  £  0)  because  Rm(x<0)  =  0  in 

equations  (26)  and  (27),  even  though,  in  general,  g^  and  gz 

[equations  (39)  and  (40)]  will  be  nonzero  upstream.  As  x  is 

increased  from  -®,  everything  remains  constant  at  its  upstream 

value  until  x  =  0_.  From  x  =  C  to  x  =  0+,  the  viscous 

subshock  occu  In  (v  ,B  ,B  )  space,  this  is  a  straight  line 

x  y  z 

parallel  to  the  v  axis,  starting  at  the  upstream  point  and 

A 

terminating  at  another  point  on  f  *  0  given  by  (vx»Byi *Bzi 
where  v* c 1 .  At  this  point,  Rm  is  no  longer  zero  (by  assumption), 
and  as  x  is  further  increased,  the  trajectory  moves  away  from 
the  subshock  point  along  an  integral  curve  which  lies  entirely  in 
the  surface  f  =  0.  As  explained  in  the  Appendix,  this  "magnetic11 
part  of  the  trajectory  must  be  on  the  subsonic  sheet  of  f  =  0 
(whereas  the  upstream  point  is  necessarily  on  the  supersoric 
sheet);  the  sonic  curve  on  f  =  0  is  the  locus  of  points  of  tan- 
gency  of  lines  parallel  to  the  vx  axis.  As  x  -  +®,  the  trajectory 
approaches  the  downstream  (nodal*)  singularity,  given  by  the 


Except  in  isolated  special  cases,  as  explained  in  §  3-3, 
when  the  downstream  subshock  point  lies  exactly  on  the  saddle- 
point  separatr ix  [corresponding  to  a  particular  value  of  E* 

(=EMAX^* 
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intersection  of  the  hyperbola  gy  =  0  and  the  (vx,8y)  plane 
curve  f  =  0.  The  part  of  the  solution  trajectory  on  f  =  0 
is  only  one  of  a  family  of  integral  curves  which  densely  cover 
this  surface.  The  Appendix  discusses  the  topology  of  these 
integral  curves  and  some  typical  shapes  of  the  surface  f  =  0 
on  which  they  lie. 

Figure  1  shows  two  projections  of  the  complete  (v  ,B  ,B  ) 

x  y  z 

structure  trajectory  for  a  given  set  of  parameters:  y  =  5/3, 

M1  =  *,  By,  -  1,  Bz]  =  1,  «*,  »  1.5,  and  E*  =  -C.30.  The 
corresponding  structure  profiles  are  shown  in  figure  2. 

Figure  3  is  a  three-dimensional  sketch  of  the  behaviour  of  the 
magnetic  field  vector  throughout  the  shock  layer,  for  this 

-JL 

typical  skew  shock.  In  this  case,  the  value  of  E~  has  been 
chosen  so  that  By  changes  sign.  For  other  (positive)  values 
of  E  ,  this  component  may  increase  or  decrease  (but  remain 

x 

positive).  [The  case  E"  =  0  corresponds  to  a  sw? tch-off  shock, 
and,  as  explained  in  the  Appendix,  this  type  of  shock  is  plane- 
polarized,  and  should  therefore  properly  be  considered  as  belong¬ 
ing  to  the  family  of  obi ique  shocks  rather  than  to  that  of  skew 
shocks ' . 3 

Another  case,  corresponding  tc  the  upstream  and  downstream 
magnetic  field  vectors  being  at  right  angles  (Byj=0,Bzj=l  ),  is 
portrayed  by  its  trajectory  proj act  ions  in  figure  4,  for 

x 

E  =  Ez  >  0.  It  is  clear  from  symmetry  considerations  that, 
for  each  solution  for  a  given  positive  value  of  Ez,  there  is  a 

Vor  switch-off  shocks  there  is  no  downstream  datum  for  (y,z)- 
orientation  of  the  coordinate  system,  which  may  thus  be  chosen 
so  that  Bzj  =  0. 


FIGURE  1.  Two  projections  of  a  structure  traject 
ory  In  (v  ,  B  ,  B  )  space  shown  In  relation  to 

t\  y  £» 

the  surface  f  «®  0,  for  the  case  B.,t  «  1,  B_i  « 


In  figure 


FIGURE  3.  Three-dimensional 
magnetic  field  vector  variati 
}ayer  for  the  same  parameters 


representation  of  the 

on  through  the  shock 
as  in  figure  1. 
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m i rror- image  solution  for  the  same  B 


and  a  corresponding 


negative  value  of  Ez  with  the  same  magnitude.  [This  is  also 
true  in  the  case  of  normal  shocks  (which  also  have  B^j  =  0); 
but  in  these  cases,  &2p0  in  addition,  and  the.  structures  are 
plane-polarized.] 


in  order  to  consider  the  full  ranges  of  the  parameters 
governing  skc-w  shock  n.zgnetie  field  angles,  the  above  discussion 
implies  that  projections  of  upstream  points  need  be  taken  only 
in  the  first  quadrant  of  the  ( )  plane,  assuming  that 

i*  0.  [  By2  may  vary  ever  a  range  of  positive  and  negative 

values,  depending  on  the  value  of  E  .]  Upstream  points  with 
projections  on  the  By  axis  correspond  to  obi ique  shocks,  with 
the  origin  representing  normal  (possibly  switch-on)  shocks. 

With  upstream  piojections  on  the  Bz  axis,  only  strictly  positive 
values  of  By2  (correspond ing  to  Ez  >  0)  need  be  considered. 

I  aM  cases,  swi tch-off  shocks,  By2  =  0,  can  be  treated  as  oblique 
shocks  by  reorientation  of  the  coordinate  system. 


In  considering  typical  examples  of  skew  shocks,  the  remainder 

of  this  paper  will  be  concerned  only  with  the  single  set  of  upstream 

angle  conditions:  Byj  =  1,  Bzj  =  1,  since,  clearly,  the  large 

2 

number  of  input  parameters  (y,  Mj,  B  j,  Bzj,  ,  and  E7) 
makes  more  complete  cataloguing  out  of  the  question. 


3  -  2  Variation  with  E 

If  different  values  of  E  are  taken,  holding  the  other 
parameters  constant  the  various  resulting  hyperbolas,  gy  =  gz  =  0, 
are  qualitatively  easily  visualized;  however,  the  shape  of  the 
curve  f  =  0  changes  cons iderabley  with  E  .  It  is  therefore 
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more  convenient  to  plot  a  new  curve  in  the  (v  ,B  )  plane  which 

x  y 

is  independent  of  E  ,  giving  the  locus  of  downstream  points 
(together  with  other  extraneous  intersections  of  g^  =  0  with 
f  =  0)  as  E  is  varied  parametrically,  for  constant  values 

JL 

of  the  other  parameters.  This  is  obtained  by  eliminating  E 
and  B  t  between  equations  (46),  (47),  and  (30),  giving 


h(v  ,B  )  =  (B  -  B  |  )(b  -  [(Y-1  )a^l-"YV.x..]By!-| 

y  y  y  y  [(y-i)«xj  +  (2-y)vx] 


[(Y-l)a^  -VVX]B^| 
t  (y-  1  +  (2-y)vx3 


(Y+  1  )(vx  -  i  )(vx  -  \  ) 

axl  C  (v_  '  )aid  +  <2-v)yx] 


=  0. 


The  intersection  of  this  oval  curve  with  the  downstream  -  singular 
-  point  hyperbola  of  equation  (46)  for  a  given  value  of  E"  gives 
the  possible  location  of  the  downstream  point  in  the  (vx,By) 
plane,  provided  certain  structural  and  thermodynamic  requirements 
are  met.  This  is  shown  in  figure  5:  Y  =  5/3,  M,  =  ®,  B  ,  =  1, 

„  1  y1 

2  « 

Bz]  =  1,  and  axj  =  1.5,  for  several  values  of  E  .  The  respective 

projections  of  the  magnetic  trajectories  are  also  shown.  It  is 

clear  from  figure  5(a)  that  there  is  a  maximum*  and  a  minimum 


See  §  3.3  for  a  case  in  which  the  maximum  E"  value  is  limited  by 
the  saddle-point-separatr ix  solution. 
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value  of  E",  in  between  which  solutions  may  exist  (shown  by 
the  heavy  curve  on  a  section  of  n  =  0),  provided  they  also  satisfy 
the  thermodynamic  requirements;  e.g.,  negative  temperature  (or 
pressure)  obviously  must  not  occur  on  a  structure  trajectory. 
Outside  of  this  range  of  E  ,  no  solutions  are  structural  1 y 
possible.  The  type  of  downstream  points  corresponding  to  E  = 

U.  X 

^MAX  or  ^MIN’  w^ere  t^*e  hyperbola  g  =  0  is  tangent  to  the 
curve  h  =  0,  have  sometimes  been  referred  to  as  (generalized) 
Chapman  -  Jouguet  points,  because  the  downstream  normal  velocity 
is  just  equal  to  the  downs tream  magneto-acoustic  speed  in  these 
cases  [see  K.unkel  and  Gross  (1962)3. 

In  order  to  get  some  idea  of  the  temperature  vai lation  in 
the  magnetic  layer,  it  is  necessary  to  visualize  the  intersection 
of  a  structure  trajectory  with  the  family  of  lamellar  surfaces 

A 

of  revolution  about  the  v  axis,  given  by  T  =  const,  from 
equation  (38)*.  In  the  inf  ini te-Mach-number  case,  this  equation 
becomes 


i  (y  +  1  )  ri  i  _2  ,„2  p2  n  _  j 

a  Xv rr  VxU  Vx  2(1xl  (Br  '  Brl  }J  “  T* 


(49) 


The  intersection  of  these  surfaces  with  the  (v  ,B  )  plane  is 

x  y 

shown  in  figure  6,  which  also  includes  the  downstream  locus 
curve,  h  =  0,  the  limiting  hyperbolas,  and  a  trajectory  projection 


t 


Note  that 


this  equation 


is  also  independent  of  E". 


FiGURE  6.  Constant-temperature  surfaces  of 
revolution  shown  by  their  intersection  with 
'vx’^y^  P*ane  (dotted  curves),  in  relation 
to  the  downstream  locus  h=0,  and  the  limit¬ 
ing  hyperbolas. 
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corresponding  to  the  E  value  (=  E^qj»  say)  for  which  T2(E  ) 
has  its  maximum  value  for  this  set  of  the  other  parameters. 

In  addition  to  section  curves  for  the  surface  f  =  0  (which 
contains  the  upstream  point)  and  for  the  surface  f  =  1  (which 
contains  the  downstream  sutshock  point),  section  curves  are 
also  given  for  the  surfaces  f  =  -1,  for  reference,  and  T  =  Tj^. 
The  latter  is  clearly  given  by  that  value  of  f  for  which  the 
section  curve  is  tangent  to  h  =  0  in  the  range  where  structures 
are  possible. 


Typically,  figure  6  shows  that  the  structurally  -  allowable 
range  on  h  =  0  lies  well  inside  the  f  =  0  section  -  curve 
parablola,  and  that  all  the  trajectories  lie  inside  the  corres¬ 
ponding  paraboloid  of  revolution.  On  a  trajectory,  T  jumps 
from  0  to  1  across  the  subshock,  and  then  varies  continuously 
towards  its  downstream  equilibrium  value  (although  not  necessarily 
monotonical 1 y ) . 


Note  that  the  parabola  for  T  =  0  intersects  the  vx  axis 
2  2 

at  vxQ  =  l+jaxjBrj,  and  the  corresponding  paraboloid  meets 

2  2  2 

the  (By,Bz)  plane  at  Bf  =  Brj  +2/axj.  For  reference,  the  point 

corresponding  to  the  limiting  (maximum)  value  of  f  for  which 

there  are  real  surfaces,  T^j^,  has  also  been  shown.  This  occurs 

on  the  vx  axis  at  vx  =  ivXQ,  giving  the  value  = 

12  2  * 

-g(Y+l )  Vxc/(Y_1)*  *n  9eneraU  Tj^ax  (which  must  occur  on  h  =  0) 

A 

will  be  less  than  T, 

L  I  M. 

A 

It  has  been  found  that  T^x  usual  Iv  occurs  at  points  having 
values  of  B  j  near  zero.  By  inspection  of  equation  (49),  it  is 


clear  that  this  corresponds  to  most  of  the  magnetic  energy  hav¬ 
ing  been  converted  into  thermal  energy  (through  Joule  dissipation). 
The  fact  that  By2  is  not  exactly  zero  for  is  due  to  the 

coupling  with  the  flow  kinetic  energy  through  the  vx  terms  in 
equation  (49). 

In  skew  shocks  in  which  Bz?  =0(1),  ai)  f2  values  are 
usually  greater  than  the  viscous  temperature  T  =  1,  because  the 
transverse  magnetic  field  magnitude  has  decreased,  thus  creat¬ 
ing  thermal  energy  downstream.  However,  in  some  cases  (especially 
for  oblique  or  almost  -  oblique  shocks  operating  near  the  upper 
Chapman  -  Jcuguet  point),  the  transverse  magnitude  may  actually 
increase,  and  T2  will  be  slightly  jess  than  the  viscous 
temperature.  Thus  in  shock  -  tube  experiments  operating  in  the 
Chapman  -  Jouguet  mode,  if  one  of  the  goals  is  the  achievement 
of  hi-h  dov/nstream  temperatures  by  making  use  of  the  "magnetic 
dumping1'  effect  of  Joule  heating,  it  is  clearly  desirable  to 
try  to  produce  the  lower  Chapman  -  Jouguet  condition  rather 
than  the  upper  one. 

3.3  Variat ion  wi th 

For  other  values  of  a*,,  the  behaviour  is  qualitatively 
similar  to  the  case  for  a*.  =  1.5,  given  abo/c,  over  a  finite 
range  of  .  As  the  magnetic  pressure  ratio  is  increased, 
the  oval  curve  of  figure  5  at  first  elongates  in  the  vx  direc¬ 
tion  and  shrinks  in  the  By  direction,  bodily  moving  towards 
larger  vx  values.  At  higher  ax1  values,  it  shrinks  in  its 


32 


vx  dimension  as  well.  For  all  skew  shocks,  there  is  an  upper 
limit  on  at  which  this  oval  curve  disappears  entirely 

(for  oblique  and  normal  shocks,  there  is  no  upper  limit).  How¬ 
ever,  wel 1  before  this  limit  is  reached,  certain  structural 
properties  limit  the  range  of  E*  values  to  lie  in  a  restricted 
range,  not  extending  to  one  or  both  of  the  tangent  points  (i.e., 
one  or  both  of  the  Chapman- Jouguet  points  are  excluded).  This 
is  because  the  topology  of  the  surface  f  =  0  is  such  that, 
for  E*  values  outside  of  the  restricted  range,  the  trajectory 
passing  through  the  subshock  po*nt,  instead  of  travelling  on 
the  subsonic  sheet  directly  to  the  nodal  singularity,  runs  into 
the  sonic  curve,  at  which  point  it  must  leave  the  surface,  fly¬ 
ing  off  to  V'x  —  [This  phenomenon  shows  up  in  the  numerical 
integrations  as  vx  becoming  imaginary  in  the  subsequent  in¬ 
crement  in  x.] 

Figure  7  shows  a  set  of  trajectories  for  ax^  =  2.3.  The 
maximum  E*  value  is  not  restricted  to  1 ?e  below  the  tangent 
{ Chapman- Jouguet )  value,  and  the  behaviour  is  similar  to  that 
of  smaller  magnetic  pressure  ratios  until  E*  reaches  a  value 
near  -0.33  (compared  v/ith  the  lower  tangent  value  of  -0.444). 
Figure  8  shows  why  lower  E*  values  are  impossible. 

Note  in  figure  7  the  occurrence  of  overall  expansion  shocks. 
After  the  sudden  gasdynemic  compression  through  the  viscous  sub¬ 
shock,  MHD  interaction  expands  the  gas,  as  magnetic  energy  is 
converted  into  directed  kinetic  energy.  The  density  profile 
and  corresponding  magnetic  field,  temperature,  and  entropy  plots 
are  shown  in  figure  3.  It  should  be  noted  that  these  expansion 


IE  7.  Downstream  locus  curve  and  trajec 
Ithin  the  allowable  range.  In  this  case 
:angent  value  because  of  structural  rest 


and  ^  JAN  *  t*ie  son*c  curve  the  entropy  is  stationary,  and 
ectories  either  leave  (entropy  maxima)  or  approach  (entropy  m 
the  surface  tangentially  along  straight  lines  parallel  to  the 
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shocks  also  develop  relatively  large  downstream  temperatures 
due  to  Joule  heating.  In  such  cases,  it  is  evident  that  the 
expansion  is  not  at  the  expense  of  thermal  energy. 

2 

At  still  larger  values  of  ax^ »  there  is  a  restriction  on 

the  maximum  E*  value  as  well,  as  shown  in  figure  10  for  a^,  = 

2 

2.9.  There  are  no  solutions  above  ax^  =  2.92. 


Structural  restrictions  of  a  different  nature  occur  at 


values  of  numerically  near  v*  (  =  0.25  in  this  infi- 

ni te-Mach-number  case).  This  concerns  the  relative  position 

of  the  separatrix  going  into  the  saddle  point  which  occurs  on 

2 


the  subsonic  sheet  of  f  =  0  for  these  smaller  values  of  a' 


xl 


In  some  cases,  the  upper  tangent  value  of  E*  is  such  that 
this  separatrix  does  not  enclose  the  subshock  point,  as  shown 

2  x 

in  figure  11(a)  for  aTj  =  0.27,  and  E^AY  is  restricted  to  a 


MAX 


smaller  value,  as  shown  in  figure  11(b).  Structures  of  this 
latter  limiting  type  are  actually  composed  of  two  shocks:  a 
skew  shock  (including  the  viscous  subshock)  leading  to  condi¬ 
tions  given  by  the  saddle-point  singularity,  followed  (at  an 
arbitrary  distance  downstream)  by  a  slow,  completely  subsonic, 
MHD  ob  1  i que  shock  from  the  saddle  point  to  the  nodal  singulari¬ 
ty.  The  isolated  saddle-point  solution  is  the  analogue  of  the 
(magnetically  trivial)  gas  shock  appearing  in  the  case  of  pure¬ 
ly  obTi^u©  ionizing  shock  structures  which  are  sub-Al fvenic 
downstream.  [Clearly,  there  is  a  unique  value  of  E*  asso¬ 


ciated  with  this  type  of  solution  —  a  common  feature  of  all 


FIGURE  10.  Structure1  restrictions 
on  E^x.  In  (a),  the  trajectory 
runs  into  the  sonic  curve  because 
the  downstream  point  is  fairly 
high.  In  (b),  a  solution  is 
possible  for  a  smaller  E  value 
due  to  the  lower  downstream  point 
and  the  fact  that  the  kidney¬ 
shaped  sonic  curve  projection 
has  not  extended  very  far. 

In  (c),  for  a  smaller  E  value, 
although  the  downstream  point 
is  sti’l  lower,  the  son:c  curve 
has  elongated  so  far  chat  a 

solution  is  not  possible.  For 

2  .  . 

larger  values  of  conditions 

(a)  and  (c)  coalesce. 
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E  11.  Structural  restrictions  on  for  axj  near  v^.  In  (a), 

ubshock  point  lies  on  an  unacceptable  trajectory.  in  (b),  L  =  E^AX; 
is  the  limit  solution  where  the  subshock  point  is  on  the  separatrix. 
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shock  waves  involving  downstream  saddle  points . ]  However, 
there  are  no  skew  shock  structures  analogous  to  the  super- 
Alfvenic  gas  shocks  (because  skew  shocks  must  always  be  sub- 
Alfvenic  downstream),  and  so  there  is  a  distinct  lower  limit 
to  axj  (upper  limit  to  the  shock  Alfven  number  Aj )  beyond 
which  skew  shocks  do  not  occur,  for  a  given  set  of  upstream 
field  angles. 

These  properties  can  be  summarized  by  making  graphs  of 
the  ranges  of  E*  and  the  downstream  quantities  By2,  VX2» 

A 

and  f2,  plotted  against  ,  or,  preferably,  against 
Aj  (=  a~j  ),  which  gives  a  direct  measure  of  shock  speed.  For 
comparison  with  earlier  work  by  Taussig  (1965,  1967)  on  nor¬ 
mal  and  oblique  shocks,  the  range  of  the  downstream  Alfven 
number  A2  =  vx2,/|ip2  ^  “  ^1*-/vx2  lS  a^so  Pl°ttec*  against 

A^ .  Figures  12(a)- (e)  show  these  ranges.  The  details  of  E*, 

*  ic-  i 

vx2,  and  T2  for  A^  near  vx~2  are  shown  on  a  finer  scale 
in  figures  13(a),  (b),  and  (c). 

In  addition  to  the  Chapman- Jouguet  and  structurally- 
limited  boundaries  in  figures  12  and  13  (solid  curves),  the 
tangent  values  are  shown  (dotted),  and  also  the  locus  of  maxi¬ 
mum  temperature  points  (light  dashed).  The  saddle-point  solu¬ 
tions  are  shown  by  the  heavy  dashed  curve.  Expansion  shocks 
occur  in  the  shaded  region.  Switch-off  shocks  have  E*  -  By2  = 

A 

The  plots  of  vx2,  A2,  a.td  T2  can  be  visualized  as  pro¬ 
jections  of  three-dimens ionul  surfaces,  the  other  coordinate 
being  E*.  In  this  way,  the  double-valued  regions  may,  perhaps. 


fi  <  r*‘ 
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FIGURE  13.  Detailed  behaviour  of  the  al lowable- range 
plots  near  Aj  =  (v*)"?.  The  dashed  curves  show  the 
locus  of  saddle-point  solutions  cf .  figure  11  (b)  . 
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be  more  readily  comprehended. 

Finally,  it  should  be  noted  that,  under  optimum  condi¬ 
tions,  significant  Joule  heating  can  occur,  giving  maximum 
downstream  temperatures  up  to  about  13  times  as  large  as 
the  viscous  subshock  temperature,  for  this  particular  set  of 
upstream  field  angles  (Byj  =  Bzj  =1).  This  is  because  mag¬ 
netic  energy  has  been  converted  into  thermal  energy,  the  effect 

2 

becoming  rapidly  stronger  as  the  magnetic  pressure  ratio  axj 
is  increased  (or  Alfven  number  decreased)  until  struc¬ 

tural  limitations  are  met. 
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Concl us  ion 


The  most  general  form  of  ionizing  shock  wave  involves  a 
magnetic  field  vector  whose  transverse  component  may  rotate 
while  changing  magnitude  throughout  the  shock  layer.  These 
are  called  skew  shocks  because  the  upstream  and  downstream 
magnetic  field  vectors  and  the  shock  wave  normal  do  not,  in 
general,  all  lie  in  a  single  plane.  This  p  lenomenon  arises 
because  the  electric  field  in  the  nonconducting  upstream  region 
has  a  certain  amount  of  freedom  in  both  magnitude  and  di rect ion. 
Much  has  been  written  previously  on  the  degree  for  lack)  of 
freedom  of  the  transverse  electric  field  magnitude  in  ionizina 
shock  waves.  There  seems  to  have  been  a  certain  amount  of 
hesitancy  in  being  prepared  tc  assign  the  electric  field  the 
status  of  a  parameter,  and  many  studies  have  apparently  been 
motivated  by  an  attempt  to  find  conditions  requiring  a  unique 
value  of  the  Ji  -  field.  The  philosophy  of  the  present  approach 
has  been  to  treat  the  E-  field  as  a  (free)  parameter,  allow¬ 
ing  the  detailed  structure  to  determine  whether  or  not  there 
are  unique  values  associated  with  certain  kinds  of  solution. 

As  is  wel 1  known  in  the  theory  of  nonlinear  two-point  boundary- 
value  problems,  parameter-un iqueness  is  always  associated  with 
the  requirement  of  a  solution  trajectory  to  lie  exactly  on  a 
saddle-point  separatrix.  Thus,  in  ionizing  shock  v/ave  theory, 
those  structures  which  involve  a  downstream  saddle  point  will 
be  associated  with  a  unique  value  of  the  transverse  E-  field; 
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i.e.,  the  field  will  be  structurally  determinate.  [The 

downs tream  singularity  is  the  determining  one  for  ionizing 

shock  waves  because  upstream,  there  is  no  singular  point  in 

phase  space  but,  rather,  a  whole  region  (lying  on  f  =  0)  in 

which  o  =  0.3  The  saddle  point  may  lie  in  the  (v  ,B  )  plane 

x  y 

as  is  the  case  for  oblique  (and  normal)  ionizing  shock  v/aves, 
or  (for  Pm  =  0)  on  the  f  =  0  surface  [as  seen,  for  example, 
by  projection  onto  the  (8^,  Bz)  plane]  in  the  case  of  skew 
shocks.  The  former  may  have  structures  for  super- Al fven ic  as 
well  as  sub-Alfvenic  saddle  points,  but  the  latter  exist  only 
for  sub-Alfvenic  saddle  points.  Transverse  shocks  always  in¬ 
volve  a  super- Al fven ic  downstream  saddle  point  in  the  (v  ,  B  ) 

x  y 

plane.  In  zero-Pm  theory,  all  (vx,By)  -  saddle-point  structures 
are  magnetically  trivial  gas  shocks. 

Much  more  interesting,  however,  is  the  range  of  structures 
possible  when  the  downstream  point  is  a  (necessarily  sub-Alfvenic) 
node.  In  this  case  the  JE- field  is  a  free  parameter,  and  can 
range  over  a  continuum  of  values  between  limits  determined  by 
other  kinds  of  topological  features  of  the  structure.  In 
zero-Pm  skew  shocks,  this  node  lies  in  the  subsonic  sheet  of 
the  f  =  0  surface.  [In  fact,  this  singularity  (which  is  cha¬ 
racterized  by  vx  <  cg,  the  slow  magneto-acoustic  speed)  is  also 

nodal  in  the  (v  ,  B  )  plane;  however,  for  zero  Pm,  the  form 

x  y 

is  degenerate,  and  all  trajectories  enter  along  the  curve  f  =  0, 
except  for  one  straight  line  parallel  to  the  vx  axis.] 
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One  difficulty  In  studying  ionizing  shock  waves  is  in 
cataloguing  the  large  number  of  different  types  of  solutions 
for  ranges  of  the  f ? ve  governing  parameters:  Mach  number, 

Alfven  number,  two  field  angles,  and  the  transverse  electric 
field.  This  analysis  has  been  restricted  to  M-j  =  ®,  and 
most  of  the  solutions  are  given  for  8^  =  B ^  =1;  thus,  Aj 

A 

and  E  (actually,  E*  =  E, +  A,  B  .)  are  the  main  parameters 

z  4  y  i 

of  interest. 

Finally,  it  should  be  stressed  that  this  analysis  has 
studied  the  magnetic  structure  of  skew  shocks.  This  is  made 
possible  by  the  coordinate  transformation  of  §  2.3,  which  de¬ 
couples  variable  ionization  and  nonequilibrium  electron  tempera¬ 
ture  effects  on  the  electrical  conductivity.  The  structure 
of  these  variables  can  now  be  studied  separately,  with  the 
known  magnetic  structure  being  used  as  an  input  (through  Joule 
heating  and  the  MID  interaction  on  the  fluid  variables).  The 
two  classes  of  problem  are  finally  combined  by  making  the  in¬ 
verse  coordinate  transformation  (a  single  straight-forward 
numerical  integration).  Work  is  proceeding  along  these  lines. 
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The  extensive  numerical  calculations  were  expertly  performed 
by  June  Clean»an  with  the  aid  of  an  electronic  digital  computer. 


Integral  curves  on  f  =  0 


The  shape  of  the  surface  f  =  0  is  defined  by  a  family 
of  circles  v  =  const.,  whose  centres  lie  on  the  centreline 

✓v 

hyperbola:  B  =  B  .  (v  ),  B  =  0.  The  intersection  of  the 
y  y<£  x  z 

surface  with  the  {vx,  B^)  plane  forms  the  curve  fQ  =  0, 
the  detailed  shape  of  which  varies  considerably  depending  on 


the  value  of  E*,  and  the  location  of  the  asymptote  of  the 

ii 


centreline  hyperbola,  vA  =  (y-I)^.  /y,  relative  to  v  =  v* 

A  A  I  A  A 


and  vx  =  1. 


There  are  so  many  topologically  different  characteristic 
shapes,  that  it  is  not  feasible  to  indicate  all  of  them  here. 
Instead,  a  few  cases  will  be  chosen  which  most  clearly  demon¬ 
strate  the  qualitative  behaviour  of  the  integral  curves  lying 
on  f  =  0. 


This  surface  can  be  considered  to  be  divided  into  two 

“sheets"  by  the  space  curve  lying  on  it  formed  by  the  locus 

of  points  of  tangency  of  lines  parallel  to  the  v  axis.  Since 

grad  f  Tin  (v  ,  B  ,  B  )  space}  must  be  perpendicular  to  the 
x  y  z 

v  axis  on  this  curve,  it  is  given  by  the  intersection  of  f  =  0 
with  the  surface 


f  af  ) 

V8Vx'vbz 


(A. 


i} 


A. 2 


from  equations  (29),  (30),  (31),  and  (33),  where  M  Is  the 
local  Mach  nunber,  given  by 


«2 


vx  vx2l 


(A. 2) 


since  the  ordinary  sound  speed  a  is  proportional  to  T  2. 

It  can  be  shows  that  M  is  a  monotonical ly  increasing  function 
of  vx,  for  a  fixed  jjlj ,  thus  the  two  sheets  of  f  =  0  can 
be  labelled  as  "supersonic"  (corresponding  to  larger  vx  values) 
and  "subsonic",  being  sparated  by  the  "sonic  curve",  on  which 
M  =  1 .  As  with  ordinary  gas  shocks,  an  acceptable  viscous  sub¬ 
shock  trajectory  must  have  its  upstream  end  or  the-  supersonic 
sheet  and  its  downstream  end  on  the  subsonic  sheet.  The  sub¬ 
sequent  magnetic  portion  of  the  trajectory  must,  if  it  is  to 
be  a  solution  trajectory,  lie  entirely  on  the  subsonic  sheet, 
because  trajectories  on  f  =  0  which  reach  the  sonic  curve 
subsequently  fly-off  to  vx  =  +«*>.  This  behaviour  is  similar 
to  that  of  siow  oblique  shocks  in  a  pre- i on i zed  gas  [especially 
in  cases  of  tensor  conductivity  (Leonard,  i9u6)3,  which  always 
have  their  magnetic  trajectories  on  the  subsonic  sheet  of 
f  =  0,  even  when  there  is  no  subshock  (Hi  ^  1).  Fast  oblique 
or  normal  shocks,  however,  have  their  magnetic  structure  on  the 
supersonic  sheet,  upstream  (where  ^  0)  of  a  suoshock,  if 
one  exists. 

Figure  A.1  represents  the  case:  y  =  'V 3 , 


FIGURE  A.l.  Representation  of  the  surface  f  =»  0  by  projection  of  section  curves 
v  =  const,  on  the  subsonic  sheet.  In  (a),  these  are  straight  lines  terminating  on 
f0  =  0  and  the  sonic  curve.  In  (b),  they  are  circular  arcs  terminating  tangent¬ 
ially  at  the  sonic  curve.  Projections  of  the  solution  trajectory  are  also  shown. 


Byj  -  1,  Bz|  =  1,  a Ji  =  0.5,  E*  =  0.25.  [For  M,  =  -  and 

Y  =  ^/3,  equation  (A. 2)  is  written  M2  =  l6v2  /  ( 5 T ) . 3  in  the 

y\ 

(vx,By)  projection  of  figure  A. 1(a),  the  curve  fQ  =  0  is 

shown,  together  with  the  downstream  hyperbola  gy  =  0  and  the 

centreline  hyperbola  B  =  B  .(v  ).  The  arrow  directions  on 

/  Vt  x 

fQ  =  0  correspond  to  the  sign  of  gy  in  the  three  different 
regions  determined  by  the  branches  of  the  downstream  hyperbola. 
Since  gz  =  0  in  the  (vx,By)  P^ane»  this  determines  those 
integral  curves  which  lie  in  this  plane.  In  figure  A. 1(b), 
the  circles  vx  =  const,  on  the  subsonic  sheet  of  f  =  0  are 

shown;  i.e.,  this  represents  the  view  as  seen  from  the  origin, 

looking  in  the  positive  v  direction. 

A 

To  find  the  shape  of  the  remaining  integral  curves  which 
lie  on  the  subsonic  sheet,  it  is  convenient  to  study  their 
direction-field  projection  in  the  (By,  Bz)  plane.  From  equa¬ 
tions  (39)  and  (40),  it  is  seen  that 


'  VX)B2 


[bv-e*/ (<,  -vxn 

B_ 


Note  from  this  equation  that  if  v 


(A. 3) 


=  const,  [a  circle  in 


figure  A. 1(b)?,  a  given  value  of  dBy  /  dB2  =  5,  say,  gives 

the  straight  line  By-E*/  (a^  -vx)  =  SBZ,  which  itself  has 

slope  S.  Where  this  straight  line  intersects  the  circle 

v  =  const.,  the  integral  curves  at  those  points  have  slope  S. 

In  other  words,  at  any  point  (v^^,  B^,  B^^)  on  the  sur- 

x  y  z 

face  f  =  0,  the  (By,  Bz)  projection  of  the  integral  curve 
passing  through  that  point  is  directed  along  a  straight  line 
joining  the  point  (By^,  B^^)  with  the  point  (By9^,  B29^), 
where  By9^  =  E*/  (a^  -  v^^),  and  Bz9^  =  0.  The  sense  of 


fhe  sense  of 


the  arrows  representing  dx  >  0  is  determined  by  djB^j  /  dx  <  0 
as  (a^.  -v^)  ^  0  [or  by  continuity,  on  the  circle  v^^  = 

X I  X  X 

a^.3,  and  by  dB  /  dx  ^  0  as  g  ^  0,  if  B  =  0  fas  with 
xi  y  y  z 

oblique  and  normal  shocks!.  To  make  this  clearer,  figure  A. 2 
shows  four  cases  depending  on  the  location  of  the  "fan  point" 
(B^9),0)  : 

(a)  inside  the  circle, 

(b)  on  the  circle, 

(c)  outside  the  circle,  and 

(d)  at  infinity  (vx^  =  ct^j). 

fThe  sense  of  the  arrows  has  been  assigned  arbitrarily.! 


Given  only  the  curve  fQ  =  0,  it  is  a  very  straight-for¬ 
ward  matter  to  make  a  graphical  construct  ion  of  the  complete 
direction  field,  from  which  the  integral  curves  can  be  deter¬ 
mined,  as  shown  in  figure  A. 3(b);  the  (v  ,  B  )  projection  of 

x  y 

figure  A. 3(a)  is  found  by  cross-plotting. 


SSHi 


Note  that  the  downstream  point  is  a  node,  whereas  the 


A  C 


other  intersection  of  gy  =  0  with  fQ  =  0  is  a  saddle 
point.  The  respective  separatrices  have  been  shown  in  slight¬ 
ly  heavier  curves. 

For  given  upstream  conditions  [in  the  present  case,  the 
point  (1,  1,  1)]  the  solution  trajectory  jumps  from  the  su¬ 
personic  sheet  to  the  subsonic  sheet  at  (v*,  1,  1),  from 
which  point  on  it  follows  the  integral  curve  passing  through 
that  point  until  it  reaches  the  downstream  singularity.  Note 
that  if  B  ^  =  0  (oblique  or  normal  shocks),  the  complete 

trajectory  lies  entirely  in  the  (v  ,  B  )  plane,  the  magnetic 

x  y 

portion  following  fQ  =  0,  which  in  this  case  passes  through 
the  point  (v*#  1,  0).  Also,  if  E*=  0,  the  surface  f  =  0 
becomes  a  surface  of  revolution,  and  all  trajectories  are 
radial  lines,  as  shown  in  figure  A. k.  Clearly,  there  is  no 
distinction  between  oblique  and  skew  shocks  in  this  case; 
i.e,,  all  switch-off  shocks  with  the  same  |JBj  j  are  plane- 
polarized,  and  have  identical  |j5|  —  structure.  The  downstream 
singularity  in  this  case  is  a  special  type  of  completely  sym¬ 
metrical  node  ('‘sink’1);  for  small  |E*|  values,  the  singulari¬ 
ty  has  the  more  standard  form  of  figure  A.3(b)  locally,  al¬ 
though,  further  away  from  the  node,  most  of  the  trajectories 
are  fairly  straight  in  (By,  Bz)  —  projection. 

2 

Finally,  figure  A. 5  shows  the  case  for  ax,  =  0.5, 

E*  =  =  0.2976.  This  is  typical  of  the  upper  Chapman- 

2 

Jouguet  solutions  for  moderate  values  of  axj. 


The  singularity 


FIGURE  A, 4.  Plane-polarized  switch-off  shock  wave  structure  for 

O 

Note  the  unusual  circular  singular  curve  at  v  ®  a  ,  (*  0.3,  i 


IGURE  A, 5.  Trajectories  for  »  ETAN*  Note  the  hybrid  singularity  which  has 

addle-point  characteristics  on  one  side,  and  nodal  characteristics  on  the  other. 


in  this  case  has  saddi e-point  behaviour  above  the  separatrix 
and  nodal  behaviour  below.  As  mentioned  in  §3.3,  the  up¬ 
stream  point  must  lie  within  the  separatrix,  as  shown,  for 
a  Chapman- Jouguet  solution  to  exist.  A  similar  situation 
occurs  for  the  usual  type  of  Ej^,  except  the  position  of 
the  singularity  is  well  below  (and  possibly  below  = 

and  the  form  of  the  singularity  is  the  inversion  of  that  shown 
in  figure  A. 5.  This  singularity  has  some  rather  interesting 
properties  in  that  it  cannot  be  linearized  locally,  and  tra¬ 
jectories  entering  on  the  nodal  side  do  so  only  extremely 
slowly  (as  if  they  were  being  "pushed  away"  by  the  outgoing 
tendencies  of  the  saddle-point  side). 
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